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Multiband lightcurves of tidal disruption events 
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ABSTRACT 

Unambiguous detection of the tidal disruption of a star would allow an assessment of the 
presence and masses of supermassive black holes in quiescent galaxies. It would also provide 
invaluable information on bulge scale stellar processes (such as two-body relaxation) via the 
rate at which stars are injected into the tidal sphere of influence of the black holes. This rate, 
in turn, is essential to predict gravitational radiation emission by compact object inspirals. The 
signature of a tidal disruption event is thought to be a fallback rate for the stellar debris onto 
the black hole that decreases as t -5 / 3 . This mass flux is often assumed to yield a luminous 
signal that decreases in time at the same rate. In this paper, we calculate the monochromatic 
lightcurves arising from such an accretion event. Differently from previous studies, we adopt 
a more realistic description of the fallback rate and of the super-Eddigton accretion physics. 
We also provide simultaneous lightcurves in optical, UV and X-rays. We show that, after a 
few months, optical and UV lightcurves scale as t~ 5 / 12 , and are thus substantially flatter than 
the i -5 / 3 behaviour, which is a prerogative of the bolometric lightcurve, only. At earlier times 
and for black hole masses < 10 7 Mq, the wind emission dominates: after reaching a peak 
of 10 41 - 10 43 erg/s at roughly a month, the lightcurve decreases steeply as ~ t 2 6 , until 
the disc contribution takes over. The X-ray band, instead, is the best place to detect the £ -5 / 3 
"smoking gun" behaviour, although it is displayed only for roughly a year, before the emission 
steepens exponentially. 
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1 INTRODUCTION 

Luminous flares from quiescent (non-AGN) galaxies are often 
interpreted as arising from the tidal disruption of stars as they 
get close to a dormant supermassive black hole in the centre 
of the galaxy dKomossa & Bade! 1 19991: iGezari et alj l2008t 120091 : 
lEsauei et al]|2008l : ICa ppelluti et alj |2009h . If this interpretation is 
correct, such events represent a very important probe for black 
holes in quiescent galaxies. In principle, detailed modeling of in- 
dividual events could yield constraints on black hole properties. 
More ambitiously, a measurement of the rate of events in galax- 
ies of different types would provide information on the stellar 
dynamics at parsec scales, that is responsib le for injecting stars 
into orbits that result in tida l disruption (e.g. JSver & Ulrnen l 19991 : 
Mago rrian & Tremainejll999T) . Knowledge of these stellar dynami- 
cal processes is of broad interest, since similar processes contribute 
to the rate of at which compact objects are capt ured by (and sub- 
sequently merge with) the central black hole ( Sigur dsson & Reesl 
1 19971) . These extreme mass ratio inspirals are important targets for 
future gravitational wave experiments. Similar processes also occur 
for lower mass objects, such as in compact binary systems, host- 
ing i) a stellar mass black hole and a white dw arf (Rosswo g et al.l 



The pionee ring works by lLacv et alj Jl982h ; iReej dl988h : 
I PhinnevI J 19891) and, fr om a numerical point of view, by 
lEvans & Kochanekl Jl989h have set the theoretical standard for the 
interpretation of such events. In particular, a distinctive feature of 
this theory is that, to a first approximation, the rate at which the 
black hole is fed with the stellar debris after the tidal disruption 
(the "fallback rate" Mft,) should decrease afl Mfb(t) oc t~ 5/3 . 
If the fallback rate can be directly translated into an accretion lu- 
minosity, one could then expect the light curve of such events to 
follow the same behaviour with time. Indeed, a t^ 5 ^ 3 light curve 
is generally fitted to th e observed luminosities of events interpreted 
as stellar disruptions jGezari et alj|2008l 120091 : lEsquei et alj|2008l : 
ICappelluti et alj|2009l) . In practice, the t" 5/3 evolution of the fall- 
back rate is only valid several months after the event, and the initial 
evolution generally shows a relatively gentle rise w ith time, de- 
pende nt on the internal structure of the disrupted star dLodato et al.1 
l2009h , 

However, it is not obvious a priori that one can make such 
a direct translation of fallback rate into observed luminosity, for 
a number of reasons. Firstly, for small black hole masses and 



2008); or ii) two white dwarfs (Rosswog 2007; Rossi & Begelman 



2009) 



1 Note that the original paper by Rees 198^ quotes a t 5 / 2 dependence, 
later corrected to t~ 5 / 3 by Phinnev 1989. 
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at optical frequencies IStrubbe & Ouataertl J2009h showed that the 
emission can contain at early times a significant component aris- 
ing from a radiatively driven wind, consequence of an inital 
super-Eddington fallback rate. Additionally, the unbound stellar 
debris might also give some contribution at optical wavelengths 
jKasen & Ramirez-Rui j I20I0I) . Secondly, observations are typi- 
cally made in some narrow bands. Since the system is intrinsi- 
cally time variable, this implies that the monochromatic lightcurves 
should not necessarily scale at all times with the bolometric curve. 

In this paper, we consider systematically both effects and we 
show that in general the expected monochromatic lightcurve in 
most bands does not scale as t~ 5 ^ 3 . When the luminosity is dom- 
inated by the accretion of the fallback material, at optical and 
UV frequencies the lightcurve is proportional to t~ 5 ^ 12 during 
an extended period of time. Only in X-rays one should expect a 
lightcurve close to the naive expectation and only for a brief period 
of time. Additionally, during the first year or so since disruption, the 
wind emission may further complicate the shape of the lightcurve 
in the optical and in the UV bands. We show that the strength of this 
contribution depends on the internal structure of the disrupted star 
and on the fraction of mass that the black hole can accrete, when 
fed at a super-Eddington rate. 

The paper is organized as follows. In section[2] we review the 
basic dynamics involved in a tidal disruption event. In section [3] 
we demonstrate analytically what are the expected monochromatic 
lightcurves of the disc and of the wind emission, separately. In sec- 
tion[4] we discuss the full lightcurve in several bands, as a function 
of the main parameters of our model. In section [5] we discuss our 
results and draw our conclusions. 



2 BASIC DYNAMICS 

Let us consider a star with mass M* = m.M© and radius R, = 
x+Rq which moves in a highly eccentric orbit under the sole influ- 
ence of a supermassive black hole with mass M = 10 6 MqMq. If 
the distance of closest approach r p to the black hole is within the 
tidal sphere of radius 



n 



i). 17 \r ['— 



1/3 



(1) 



then the star self-gravity is not able to counteract the tidal pull and 
the star is torn apart. We define the penetration factor as the ratio of 
the two distances: 



(2) 



When needed, our standard choice of parameters for the stellar 
properties is ft — 2* = m, = 1. 

For simplicity, we will assume the star trajectory to be 
parabolic. However, the dynamics described here is more gener- 
ally valid for stars whose orbital energy is much smaller than the 
energy spread caused by tidal forces. The stars that enter the tidal 
sphere are generally expected to satisfy this requirement, since they 
should have very eccentric orbits. 

At distances r larger than r t , R*/r < (M*/M) 1/3 « 1CT 2 
and the star behaves effectively as a point mass in the black hole 
gravitational field. However, when the star approaches its pericenter 
r P w rt, the star size becomes important. The different distances to 
the black hole at which different fluid elements lie cause a sizeable 
spread Ae in specific energy e within the star, 

GM 



dLodato et ai]|2009l) . which can greatly exceed the star binding en- 
ergy Ae/(GM*/R*) = /3 2 (M/M») 1/3 . As a result of this spread- 
ing, half of the star debris remains bound (e < 0) to the black hole, 
while the rest (e > 0) is lost from the systemQ The material which 
remains bound is ejected in highly eccentric orbits and after a time 



2-kGM 
(-2Ae) 3 / 2 

21/ 2 U, 
41 Mg^mr 1 /? 



3/2 
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GM 
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(4) 



it starts coming back to pericenter at a rate Mfb, which depends 
on the distribution of specific energies. If the distribution i n energy 
dM « / de is flat, then we find the classic result (Rees 1988; Phinnev 
1989) 



Ma = M p ( 



-5/3 



The peak fallback rate 

M p = ii^ « 1.9 10 26 M- 1/2 mt f3 3 x; 3/2 g/s, 

can be compared to the Eddington accretion rate 

M Edd = 1.3 10 24 M 6 (ij/0.1) -1 g s"\ 



(5) 



(6) 



(7) 



where 77 ~ 0.1 is the radiative efficiency of the accretion process. 
We thus see that for a lO 6 A/0 black hole, the fallback rate can be 
as large as 100 times the Eddington rate, and is therefore expected 
to produce a luminous flare, which can be in principle detectable. 
Note that the ratio M p /M Bdd oc M~ 3/2 , so that the fallback rate 
is expected to be only marginally super-Eddington for a 1O 7 M0 
black hole. Although also in this case a luminous flare is expected 
to be produced, its appearence will be significantly different, as 
discussed below in Section 4. After a time 



t E dd = 760 (rj/0.1) 3/5 M- 2/5 ml /5 p- 6/5 x't /5 d, 
the fallback rate ceases to be super-Eddington. 



(8) 



This simple picture has been extensively tested numerically 



Ae : 



-R„ 



(3) 



(Nolthenius & Katz 1982; Evans & Kochanek 1989; Laguna et all 
1 1993l;lAval et alj2000l) . In particular. lLodato et all j2009h show two 
important results, that we will use in the following. First, even when 
ft = 1 and the the spin up of the star is expected to approach the 
break up velocity, eq. <[3j is a good estimate for the total energy 
spread. Second, the energy distribution dM,/de within the bound 
matter depends on the internal structure of the star (parameterized 
by its poly tropic index 7). In particular they find that dM*/de is 
constant only for an idealised incompressible star (7 2> 1). In this 
idealised case Mfb follows eq. (O at all times. In general, however, 
the resulting Mfb departs from the t ~ 5// ' 3 power-law behaviour, dur- 
ing the first few months (see their figure 10, left panel), resulting in 
a smaller peak fallback rate. The departure is more pronounced for 
more compressible stellar models. 

As matter falls back to r p , we assume that it undergoes shocks 
that convert all its kinetic energy into internal energy. Initially, 
when the accretion rate is highly super Eddington, only a frac- 
tion (1 — /out) of the fallback material forms a disc and can ac- 
crete all the way down to the black hole. The remaining part of the 



2 If the stellar orbit is not exactly parabolic, then the ratio of bound versus 
unbound material may vary. 
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stellar debris will instead leave the system under a strong radia- 
tive pressure. It is not possible to quantify analytically the parame- 
ter /out, because the amount of mass that is ejected will depend 
upon the detailed re distribution of energy and an gular momen- 
tum within the flow dBlandford & Begelmanll2004l) . Rather, sim- 
ulations are required to a ssess the rate at which a ho le can be fed in 
such extreme conditions. IStrubbe & Ouataertl d2009h assume a con- 
stant / ou t = 0.1. For moderate super-E ddington accretio n rates , 
however, numerical slim disc models by Dotan & S havivl d201Ct) 
indicate that / ou t is a growing function of A/fb/MEdd, reaching 
/out ~ 0.7 for Mfb/MEdd = 20. For higher accretion rates, they 
could not find steady state solutions. This does not mean that so- 
lutions with a different accretion geometry do not exist. However, 
it does mean that any extrapolation of their results to more vigor- 
ous accretion regimes, that may be attained as a consequence of a 
tidal disruption event, should be considered as indicative. Mindful 
of these uncertainties, we will explore the consequences of differ- 
ent prescriptions for / ou t ■ 

The material which is not ejected in the wind accretes onto 
the black hole through a relatively thick disc. Initially, the disc vis- 
cous time is much shorter than the timescale at which the material 
falls back to r p . In this phase the rate at which the black hole ac- 
cretes mass is thus the same as the rate at which the disc is fed 
from the fallback material. At ve ry late times (at t > 100 years 
after the event, e.g. IUlmeJll999t) the viscous time becomes long 
enough that the accretion rate is controlled by viscous processes 
rather than by the disc feeding. ICannizzo et al.l jl990h have dis- 
cussed these late stages and found that the accretion rate onto the 
hole scales as In this paper we will not consider these late 

stages, focussing instead on the early evolution of the system. 



3 SCALINGS FOR THE MULTIWAVELENGTH 
EMISSION 

In this section, we discuss the relevant dynamics and radiative prop- 
erties of the disc and the wind separately. In particular, we will 
focus on the temporal behaviour of monochromatic lightcurves, 
which are derived analytically assuming a pure t~ 5 ^ 3 power-law 
for the accretion rate. Therefore these scalings should be regarded 
as being only suggestive at early times (< half a year). However, 
they certainly provide a good description at times larger than one 
year after disruption, when the observed emission is dominated by 
the disc at all frequencies (see Section 4). 



3.1 Disc contribution 

In this section, we focus on the disc emission, which is expected to 
dominate at high frequencies and/or at late times. For simplicity, we 
adopt here a thin disc geometry, and defer to section|4]a discussion 
of a probably more realistic slim disc model. For a geometrically 
thin and optically thick disc, the effective temperature at a given 
radius R is related to the mass accretion rate through: 



4 3GMM 
aSBTd = ~~8ttRF~ 1 ' - 




(9) 



where M is the black hole mass, M = (1 — fout)M{h is the accre- 
tion rate, Ri n ~ 3Rs (where R B — 2GM/c 2 is the Schwarzshild 
radius of the black hole) and <tsb is the Stefan-Boltzmann constant. 
As time increases, the accretion rate decreases as t~ 5 ^ 3 . Thus, 



while the bolometric luminosity scales as M, the disc temperature 
decreases as M 1 ^ 4 . 

For a star of mass M* = Mq and radius i?* = Rg , disrupted 
by a 1O 6 M0 black hole, the tidal radius is approximately at 20i?s 
from the hole (eq.Q}. Since the disc outer radius is expected to be 
close to the circularization radius Rout = 2r p = 2r t / /3, the radial 
extent of the disc is rather small. In general, we have 



Rout 
Rin 



15 / 3- 1 M~ 2/E W 1/3 3 



(10) 



The spectrum of an optically thick disc is approximately given 
by a superposition of blackbody spectra with temperature ranging 
from Tout = T(-Rout) to T in = T(R in ). Since the radial extent of 
the disc is small, the range of disc temperatures is correspondingly 
small: 



Tin 
Tout 



§2*) ^7.6/3" 3; 

Rin } 



M, 



-1/2 



-1/4 3/4 
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If the event is observed at a given wavelength u bs, we may 
define a corresponding temperature through Wien's displacement 
law: 



To: 



1.7 lOVstf, 



(12) 



where we adopt the notation: ^ x = v bs/W x Hz. The shape of the 
lightcurve is determined by the relative size of T b s , Tin and T ou t. 
We can identify three phases: 

(a) T bs < Tout. This corresponds to the earliest times (when the 
disc is hotter), for relatively small observed frequencies. In this case 
the observed frequency sits in the Rayleigh- Jeans tail of the disc 
spectrum and thus the luminosity scales with time in the following 
way: 

T rri JlVl/4 ,-5/12 

vL v oc T d oc M ' oct ' . 



(13) 



(b) Tout < Tobs < T n . In this phase the emission at i/ bs is 
dominated by the disc radii for which T(R) ~ T Q b s , whose flux is 
close to the bolometric flux. Since the emitting area does not vary 
much with time (at most by a factor i? ut / Rin ~ 10), it is in this 
phase that we can expect a monochromatic luminosity proportional 
to t~ 5 ^. 

(c) Tin < Tobs - This phase occurs at late times, when the disc has 
cooled significantly. In this case, the observed frequency sits in the 
Wien part of the spectrum and we should expect an exponentially 
declining monochromatic lightcurve. 

The transition between the different phases occurs at times t a b and 
the, such that Tout(iab) = T b s and T in (tbc) = T b s . These two 
transition times are given by: 



tab « 6.3 10 4 /? 3/5 M 6 1/5 m^ 5 x- 6 / 5 ^- 12/5 d, 
t hc » 8.4 10V 6/6 M 6 - W 5 4 /5 WK 12/5 d- 



(14) 
(15) 



Typically these transitions occur at times much larger than tedd, 
and we have thus assumed / ou t = in Eqs. d 14fc and d 1 5 b - We 
plot in Figure Q] the relevant timescales as a function of the black 
hole mass Me for three choices of the observed frequency ^ bs = 
6.3 10 14 Hz (optical g-band), u ohs = 1.9 10 15 Hz (the GALEX 
FUV band) and u ohs = 4.8 10 16 Hz (soft X-rays at 0.2 keV). The 
solid line shows the time t m m, which marks the beginning of the 
tidal disruption event. The short dashed line indicates tab and the 
long dashed line indicates t DC . We thus clearly see that in the optical 
and in the FUV, the lightcurve initially follows a t -5//12 lightcurve 
for an extended period of time, lasting several hundred years for 
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Figure 1. Typical timescales as a function of black hole mass. The solid line shows t m i n (Eq.[4}, the short dashed line shows and the long dashed line 
shows tb c . The three panels refer to an observed frequency in the optical g-band (left), in the GALEX FUV band (center) and in the X-rays, at 0.2 keV (right). 
In the optical and FUV, for an extended period of time after the event the lightcurve is expected to scale as t~ 5 / 12 before steepening into t~ 5 / 3 . It is only in 
X-rays that the initial lightcurve follows the standard i -5 / 3 law. 



the g-band and several decades for the GALEX FUV band. Only 
in the X-rays is the initial lightcurve expected to show the standard 
£-5/3 £i ec ]i ne f or a f ew years, before turning into an exponential 
decline. Note also that the time span during which we expect a 
monochromatic t _5//3 decline gets smaller for higher black hole 
masses, at all frequencies. 

An example of a disc lightcurve at different frequencies is 
shown in Fig. [2] which corresponds to a solar type star disrupted 
by a 10 6 Mq black hole. The fallback rate is the one predicted by 
iLodato et all fe009l) for a 7 = 1.4 polytropic model for the star. 
The two red lines trace, for comparison, a t~ 5 ^ 3 and a t~ 5 ^ 12 de- 
cline. This figure confirms the simple estimates described above: 
the optical and FUV lightcurve show initially a £~ 5//12 decline, 
while the X-ray lightcurve follows the bolometric t~ 5 ' 3 curve; af- 
ter a few years the X-ray lightcurve steepens into an exponential 
decline. 



3.2 Wind contribution 

In this section, we derive the monochromatic luminosity scalings 
for the wind contribution to the emission that may arise from a 
tidally disrupted star. 

The main properties of a radiation-driven wind powered by 
fallback matter have been calculated b y Rossi & Begelman (2009 ) 
and applied to the problem at hand by Strubbe & Ouat aert (2009). 
Here we state some of their main results without derivation and we 
refer the reader to that literature for more details. In addition, we de- 
rive scalings to better illustrate the behaviour of a monochromatic 
lightcurve. 

The main assumptions are the following. From the launching 
radius tl ~ Rout, the wind expands adiabatically and radially, 
with a constant velocity « w , which we parametrize as a factor f v ^ 
1 times the escape velocity y/2GM/r\ J , 

w w « 4.4 10 9 /v/3 1/2 M 6 1/3 ^ 1/2 my 6 cm/s. (16) 

This is justifiable because R ou t is close to the sonic radius, after 
which the velocity varies only by a factor of a few. We note that 
the two wind parameters, f v and /out, should be determined and 
related by global conservation of energy and mass for the whole 
system. However, given the level of simplification of our model, we 




10 39 1 1 — 1 — ' — 

10 100 1000 10 4 

t(days) 

Figure 2. Lightcurves for the disc emission from the disruption of a solar 
type star by a 10 Mq black hole at /3 = 1. The fallback rate corresponds 
to a 7 = 1.4 polytropic model for the star jLodato et al.ll200^) . The time 
t = corresponds to the pericenter passage of the disrupted star. The solid 
line shows the bolometric luminosity, the dot-dashed line shows the lumi- 
nosity at 0.2 keV, the long-dashed line shows the luminosity in the FUV, at 
v = 1.9 10 15 Hz, while the short-dashed line shows the luminosity in the 
optical at v = 6.3 10 14 Hz. The two red lines mark the simple power laws 
expected for the bolometric luminosity (oc i~ 5 / 3 ) and for the monochro- 
matic luminosity in the optical/UV (oc t~ 5 / 12 ). 



do not attempt to calculate this relationship, which depends on the 
details of how energy and mass are actually ridistributed between 
the disc and the wind. For this reason, in the following, we will 
treat these parameters as independent. 

The radiation is advected with the flow up to the trapping ra- 
dius, where the optical depth is r ~ c/v w . Afterwards, the radiative 
transfer takes over, until photons reach the photospheric radius r p h. 
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Since c/v w m 7, the trapping radius and the photosphere are close 
enough that we may assume that the radiation keeps cooling along 
the same adiabatic curve for the whole evolution. As we show in the 
following the temperatures in the flow are high enough to justify the 
use a constant Thomson opacity for a fully ionized gas n s — 0.4 g 
cm -2 . 

The radiation temperature at the base of the wind, TL = 
T(tl) can be derived from energy conservation in the wind. If all 
the internal energy at the base of the wind is converted into kinetic 
energy at large distances the conservation law can be written as 



47rr 2 « w f ^aT 4 j ~ ^Mfb/out«w, 



(17) 



from which we get 

T L (t min ) « 3.5 10 5 (/out/v) 1/ V 1/8 M 6 - S/2 V 1 mi 7 '' 24 K.(18) 

However, the wind is highly opaque and photons are released only 
at a much larger radius, the photospheric radius 

Ij- \ ^ -i A 1 nl5 J° u * a5/2 j r— 5/6 —1 11/6 ,, n \ 

r-ph(tmm) ~ 1.4 10 —f—p M 6 ' m» ' cm, (19) 

which is ~ 100 times r^. The corresponding photospheric temper- 
ature is T ph = T L (r ph /r L )- 2/3 (/ out //v) 1/3 , given by 



1 ph 



lmmJ /? 23 / 24 * {font) 



7 \ 1/12 / M *l \ I/ 72 



K. (20) 



This gives an initial flash, most luminous in the optical band. As 
the accretion rate decreases, the photosphere sinks inwards as 



Tph — ^ph(^min) 



t 



-5/3 



and the corresponding temperature increases in time as 

25/36 

Tph — Tph(imin) I ~ 

\ t-min 



(21) 



(22) 



If initially (at t m i a ) the observed frequency v ^ s lies much below 
the black body peak frequency, then the monochromatic lightcurve 
decreases quite steeply as 

r 2 rri ,-95/36 

v obs L„ cx r ph 7 ph oc t . 



(23) 



The peak in the lightcurve is thus at t = i m i n , 

"/ f23 \ 1/12 / 311 s 1/72 

3.7 10 41 ^ 3 4 /3 97/24 f^f « ' — i 



^J£^bs, rr, /. V /. N2 

o "vb-^ph (t"nxin) ' ph (t'min) 



MJ 9 



(24) 



erg/s, 



which, in the optical band, can be 10" 2 - 10~ 3 of the Eddington 
luminosity. 

When, instead, the observed frequency initially lies in the 
Wien part of the spectrum, then the lightcurve has an initial ex- 
ponential rise. The time of the peak is reached when T p h (eg I22t 
equals T b s (Eq[T2]|. 



1.2 10 



, 3 /^v /2 7^ 1/25 ^ 51x 1/50 



v^6 8 ; v/v 2 v v /? 81 



d (25) 



Afterwards, the lightcurve decreases as eq. J23l >. The peak luminos- 
ity in this case is 

^ob s £p,i/! tmax = 47rr p i 1 (f m ax) 2 o"T' 4 bs = 

i r 1n 43 -4/5 2/15 / P \ 2/,J / /out \ ^ ,. r 16/15 , (^6) 

= 1.5 10 J ^ 16 mj [ — ) ( ysr ' C1 - /s 



If the observed frequency is in X-rays, at 0.2 keV, the peak emission 
is « 10 41 erg/s at t max « 9000 days (for / out = 0.1), which 
means that the wind contribution in the X-ray is much lower than 
the disc one (see Fig.[2]> for the whole duration of the wind t ^ t c dd- 
Note that for the example above £ max ~ 9000 days 3> t c dd so the 
peak in the soft X-ray lightcurve is never reached, because the wind 
ceases before. 

In reality, the fallbak rate may not be a pure power-law, but 
rather an exponentially increasing function, followed by a plateau 
which will then lead to the t~ 5 ^ 3 decreasing slope. This has three 
main consequences. First, the duration of the super-Eddington 
phase (thus of the wind) may be shorter. Then, lower accretion rates 
are attained with respect to the power-law extrapolation at early 
times. Finally, there may be two spectral breaks in the lightcurve, 
since at early times, when M increases, the photospheric tempera- 
ture decreases, while at late times, when M decreases, the opposite 
occurs. Therefore, it may be possible that the peak of the black body 
spectrum crosses twice the observed frequency. 



4 TOTAL LIGHTCURVE 

Having discussed the fundamental scalings of the two components 
of the emission separately, we now turn our attention to the total 
lightcurve. In particular, we discuss the consequences of relaxing 
two of our previous main assumptions: (a) that the fallback rate 
follows a simple power-law decline at all times and (b) that the gas 
fraction ejected through a radiative wind / ou t is constant, indepen- 
dent on the fallback rate. Also, we will adopt a slim disc configu- 
ration during the super-Eddington phase, where the characteristic 
disc temperature is modified as: 



OSBTd = 



SGMMf 
8itR s 



2 I 4 V»?M Edd 



,(27) 



dStrubbe & Ouataertll2009l) . where / = 1 - y/BsT/R. 

The accre- 
tion rate in the disc is given by M = (1 — /out)Mfb, where the 
fallback rate is computed either from t he analytical expre ssion (Eq. 
l[5}), or from the numerical results of Lodat oet alj d2009h . depend- 
ing on the models. From Eq. d27t we directly compute the disc 
effective temperature which allows us to compute the monochro- 
matic luminosities as a function of time. 



4.1 Modeling the fallback rate: wind suppression 

As mentioned above, eqs. (5) and l|6), only hold under the assump- 
tion that the specific energ y distribution withi n the disrupted star 
is constant. As shown by lLodato et "all j2009h . the actual energy 
distribution depends on the stellar structure and only approaches a 
constant value in the limit of a completely incompressible star. If 
the star is compressible, the net effect is a more gentle approach to 
the peak, with significant emission starting at times slightly earlier 
than t m in. The peak fallback rate is consequently decreased with 
respect to the prediction of Eq. ((SJ. For example, if the star is mod- 
eled as a polytropic sphere of index 7 = 1.4, the peak accretion rate 
for a black hole mass M = 10 6 M© is only M p w 40M Ed d, to be 
compared with M p > 100A/Edd for the case of a pure power-law 
decline. We thus expect that such a more detailed modeling of the 
fallback rate would result in a suppression of the wind contribution. 

We computed the total lightcurves in several monochromatic 
bands for the two extreme cases of an incompressible star (pure 
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Figure 3. G-band (left) and FUV (right) lightcurves of tidal disruption events, for a solar type star, fj = 1, a black hole mass M = 10 6 M Q , / ou t = 0.1 and 
/ v = 1. The black curves refer to a pure power law decline of the fallback rate, while the blue lines adopt the detailed modeling of the stellar structure of 
Lodato et al. (2009), for 7 = 1.4. The solid lines are the total lightcurves, the short-dashed line indicates the wind component, and the long-dashed line the 
disc component. In all cases, after ss 1 year the disc dominates the emission. When a more detailed model for the fallback rate is used, which takes into account 
the internal structure of the disrupted star, the wind contribution is severlely suppressed and only provides a modest contribution in the optical. Lightcurves for 
the case of a 10 5 M Q and of a 10 7 M Q black holes are shown in Figs.|4]and[5] respectively. 
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Figure 4. Same as Fig.O but for M = 10 5 Mq. 



power law decline in Mf b ) and for the most compressible star con- 
sidered by lLodato et al.l ( I2009I) . with 7 = 1.4. In particular, we 
computed the fallback rate based on their equations (6)-(9), includ- 
ing the homologous expansion by a factor ~ 2.5, as observed in 
their numerical simulations. 

Fig.[3]shows the total lightcurves (solid lines), and the separate 



disc (long-dashed lines) and wind (short-dashed lines) components 
for M = 1 6 M fl . We also adopt ,f ou t = 0.1 and f v = l for 
the wind (cf. IStrubbe & Ouataertll 20091) . The left panel shows the 
g-band lightcurve, while the right panel shown the FUV one. The 
black lines correspond to a pure powe r law decline in the fallback 
rate, while the blue lines refer to the lLodato et al. I J2009h model. 



Tidal disruption events 7 




10+° r 



100 1000 

t(days) 



io 4 



Figure 5. G-band (short-dashed line), FUV (long-dashed line) and X-ray in 
the 0.2 keV band (dot-dashed line) lightcurves, for / ou t = 0.1, f v = 1, 
a black hole mass of 10 7 M o , and our standard parameters = m„ = 
(3 = 1. The red line marks a f -5 / 3 decline. 




M/M Edd 

Figure 6. The solid line shows the relation between / ou t and Mfb / Mgdd 
that we have used to model the super-Eddington accretion with a time- 
depende nt fraction of ejected m aterial (Eq. 128) . The squares indicate the 
results of lDotan & Shavivl 1201$) . 



In all cases, after « 1 year the disc dominates the emission in all 
bands, and the total lightcurve follows the t~ 5 ^ 12 decline discussed 
in section[3] In the FUV band the disc contribution dominates at all 
times. In the g-band, instead, the wind produces a luminous flare 
when the fallback rate is modeled as a pure power-law. Note that 
ou r predicted lightcurve is s lightly different that the one presented 
in IStrubbe & Ouataer 3 J2009h (their fig. 3, middle panel). This is 
becasue they assume a wider distribution of specific energies (by 
a factor 3), and because they have artificially softened the sharp 
rise to the peak at t — t m i n predicted in this case (Strubbe, private 
communication). At late times, the wind emission declines steeply 
with time as £ -95 / 36 ; as predicted by the simple scalings of section 

El 

As expected, the effect of a more gentle approach to the peak 
fallback rate is to severely suppress the wind contribution, which, 
for a polytropic star of index 7 = 1.4, provides only a modest en- 
hancement of the luminosity in the optical, while the disc emission 
dominates at almost all times. 

We have also computed the same lightcurves for a smaller 
black hole, with a mass of 10° M©, and keeping all the other pa- 
rameters fixed. The results are shown in Fig. [4] For such a smaller 
black hole, the fallback rate is more strongly super-Eddington and 
thu s the wind contribut ion is more prominent, even when we adopt 
fhe lLodatoetai] (2009) model. In this case, the wind lightcurve in 
the pure power-law model also shows an initial rise to the peak, 
before gradually declining as i -95 / 36 . This is because T b s ini- 
tially lies in the Wien part of the blackbody spectrum of the wind 
photosphere. As the wind photosphere retreats and gets hotter, the 
observed frequency crosses the blackbody peak and moves to the 
Rayleigh- Jeans side, thus producing the peak in the monochromatic 
lightcurve (see section[3]>. 

Finally, we computed the lightcurves for a more massive black 
hole, with M = 1O 7 M0. In this case, the peak fallback rate is 



only marginally super-Eddington. Even adopting a pure power-law 
decline, w e find that the wind provi des a sma ll contribution in the 
optical (c f lStrubbe & Ouataertl2009l) . With the lLodato etafl J2009t) 
model the peak fallback rate is M p « 1.3 MEdd- The wind is 
thus further suppressed and it produces a negligible contribution 
at all wavelengths and at all times. The g-band, FUV and X-ray 
lightcurves are shown in Fig. [5] Even if here the disc dominates at 
all times, the lightcurve is not proportional to t~ 5/l3 (the red line 
in the figure) in any single band. Indeed, while in the optical and 
in the FUV the lightcurve follows t~ 5 ^ 12 , in the X-rays it initially 
grows, following the growth of M, then it drops as t -5 / 3 for a short 
period of time, before eventually declining exponentially. 

We thus see that the shape of the wind lightcurve, which is 
strongly dependent on the evolution of the fallback rate at early 
times, is in turn related to the structure of the disrupted star. More 
incompressible stars give rise to an initially larger fallback rate and 
are expected to produce a stronger wind contribution. 

4.2 Modeling super-Eddington accretion: wind enhancement 

The results presented in the previous section are in fact conserva- 
tive estimates of the wind contribution. Indeed, we have assumed 
that only 10 percent of the infalling stellar debris are ejected in 
the radiatively driven wind, independently of the fallback rate. In 
reality, the strength of the wind is expected to scale with the fall- 
back rate, in such a way that very super-Eddington inflows should 
results in a much stronger wind. As mentioned above, it is not 
easy to quantify the fr action / out . We consider here the results of 
iDotan & Shavivl (2010). who have constructed self-consistent mod- 
els for super-Eddington slim discs. These authors find that the frac- 
tion /out is a growing function of the ratio Ma /MEdd, reaching 
values of the order of / ou t ~ 0.7 for Mfb/MEdd — 20. For larger 
values of the infall rate, they do not find any steady state solution, 
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Figure 7. Optical (g-band, short-dashed lines), FUV (in the GALEX FUV band, long-dashed lines), and X-ray (at 0.2 keV, dot-dashed lines ) lightcurves 
obtain ed adopting a time-dependent model for / ou t (eq. l28h for M = 10 6 Mq (left panel) and M = 10 s Mq (right panel). We also assume the Lod ato et alj 
1 2009) fallback rate for a 7 = 1.4 polytropic star. The other parameters are = m„ = ft = 1. The red line in the left panel shows, for comparison, a t~ 5 / 3 
decline. In this case, the wind contribution is strongly enhanced and now stands out both in the optical and in the FUV. On the contrary, the X-ray lightcurve 
is at all times dominated by the disc. The slight "bends" in the X-ray lightcurves are a consequence of the fact that, during the super-Eddington phase, only a 
(time-dependent) fraction of the fallback material makes it to the hole. 



and it is thus not easy to extrapolate their results beyond that point. 
We note however , that, when the fallback rate is modeled based on 
the ILodatoe t al. (2009]) results, the peak rate is not much larger 
than 20 MEdd, and does not require an extrapolation much beyond 

their maximum value . 

iDotan & Shavivl {2010) construct four different models of 
super-Eddington slim discs, with Mfb/A/Edd = 1, 5 10 and 20, 
respectively. We have approximated these results using the follow- 
ing relation between / out and Mfb /MEdd : 



/out = — arctan 

7T 



7.5 V Mi 



(28) 



which has the required properties of approaching zero when 
Mfb = MBdd and of approaching unity when Mfb 3> MEdd- 
A comparison between the relation predicted by Eq. d28t and the 
IDotan & Shavivl ( feoioh results is shown in Fig.[§] We have experi- 
mented with several different functional relations between / out and 
Mfb / M Edd and found that, as l ong as the relation reasonably fol- 
lows the lDotan & Shavivl faOlOh data, our resulting lightcurves are 
essentially unchanged. 

Fig. 13 shows th e mono chromatic lightcurves obtained assum- 
ing the lLodato et alj |2009) fallback rate for aj = 1.4 polytropic 
star, for M = 1O 6 M (left panel) and M = 1O 5 M (right panel). 
Clearly, the wind contribution is strongly enhanced and now stands 
out at early times both in the optical and in the FUV for both 
choices of the black hole mass. On the other hand, in the X-rays, the 
disc contribution dominates at all times. After 1-3 years (depend- 
ing on the black hole mass) the fallback rate drops to sub-Eddington 
values and the optical and FUV lightcurves approach the predicted 
£-5/12 (jggijj^ w hji e m e X-rays decline exponentially. During the 
super-Eddington phase the shape of the lightcurve depends on the 



combination of various factors: the time-dependence of both Mfb 
and of /out, the spectral changes occuring as the wind temperature 
increases, and the varying relative contribution of the disc and of 
the wind to the total emission. Coincidentally, the FUV lightcurve 
for a 10 6 Mq black hole does initially show a t~ 5 ^' A decline (shown 
with a red line in Fig.[6jl before flattening into t~ 5 ^ 12 . 

Finally, we considered the case of a more massive black hole, 
with 10 7 Mq. In this case the fallback rate is never much larger 
than the Eddington rate, and thus even when using Eq. J28t for / ou t 
there is no strong wind, and the emission is at all times completely 
dominated by the disc. 



5 CONCLUSIONS 

The 'signature' of a tidal disruption event is traditionally consid- 
ered to be a characteristic evolution of the fallback rate, which 
should follow approximately a t~ 5 ^ 3 decline. Indeed, in most cases 
candidate tidal disruption events are identified by observing a t ~ 5 / 3 
decay in the luminosity ove r the course of a few mon ths at various 
wavelengths from t he UV jGezari et alj|2008ll2009l) to the X-rays 
dEsquei et a l. 2008; Cappell uti et alj|2009t) . In this paper, we have 
critically reviewed the assumptions that i) the fallback/accretion 
rate should follow at all times a i -5 / 3 decay and ii) that its be- 
haviour directly translates into the observed lightcurve decay. 

As to the former, the fallback r ate behaviour is exp ected to de- 
pend on the star internal structure ILodatoet al. 2009). Generally, 
it displays initially a gentle rise to the peak and approaches a t~ 3 ^ 3 
decay only several months after the event. In this early stages, the 
rate is super-Eddington for a black hole of M < 1O 7 M0, and the 
fraction of mass that accretes onto the hole or is lost in a radiatively 
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driven wind is also time dependent <Dotan&Shavrvll201Ch . Un- 
der these conditions, the optical and UV lightcurves may display 
an intial bump given by the wind emission, which decay rapidly 
as t~ 95 ' 36 . For more massive black holes, the wind is suppressed, 
since all or nearly all the fallback mass can be accreted, and the 
lightcurve at all frequencies is dominated by the disc emission. The 
detailed evolution and the strength of the wind emission is strongly 
dependent on the evolution of fallback rate at early times, and is 
thus very sensitive to the internal structure of the disrupted star 
dLodato et al]|2009h 

Likewise, at higher frequencies and at later times, the wind 
emission is negligible for all black hole masses. Also in this cases, 
while the bolometric curve may decline as t -5//3 , the monochro- 
matic lightcurves in general do not. On a timescale of months to 
years, the optical and UV lightcurves are expected to decrease 
much slower, as t~ 5 ^ 12 . Only in X-rays does the lightcurve ini- 
tially follow more closely the bolometric luminosity, hence show- 
ing an approximately t~ 5 ^ 3 decline, before steepening into an ex- 
ponential drop after a few years. The predicted X-ray luminosities 
in the first year after the disruption are comparable to those of bright 
quasars. Detectability is thus not an issue in the local universe, and 
a standard snapshot survey with currently flying detectors would 
give enough accuracy to constrain the slope with two epochs of 
observation. Coincidentally, for some combinations of the model 
parameters, the FUV lightcurve might exhibit a decline seemingly 
close to t~ 5 ^ 3 , although clearly this does not reflect the time evo- 
lution of the fallback rate. 

Typically, observed tidal disruption events candidates have 
UV/optical luminosities in excess of 10 43 ergs/sec and lightcurves 
which decrease more steeply that t~ 5 ^ 12 . This indicates that even 
in the FUV the wind contribution might dominate at early times, 
consistent with our models which adopt a time dependent ejected 
mass fraction (see, for example, Fig. 7, left panel). We conclude 
that a more sophisticated — ideally multi-wavelength — modelling 
of the data than done so far is required to pin down a tidal disrup- 
tion event and derive the system parameters. The X-ray band seems 
best suited to identify a tidal disruption event through the fallback 
decay rate, while early optical and UV observations can catch at 
early times the wind emission, which would also teach us about 
the ill-known super-Eddington accretion physics and possibly even 
about the structure of the disrupted star. 
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